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We develop our previous study of the transition to deconfined quark phase in neutron stars, including the 
interaction in the quark equation of state to the leading order in the perturbative expansion within the confine¬ 
ment density-dependent mass model. Using the Gibbs conditions the hadron-quark mixed phase is constructed 
matching the latter with the hadron equation of state derived from the microscopic Brueckner-Hartree-Fock ap¬ 
proximation. The influence of quark interaction parameters on threshold properties and phase diagram of dense 
neutron star matter are discussed in detail. We find that the leading-order quark interaction expands the density 
range of the mixed phase, pushing forward the disappearance of the hadron phase. Moreover, since the equa¬ 
tion of state could turn out to be stiffer, a high-mass hybrid star is possible with mixed-phase core with typical 
parameter sets. 
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I. INTRODUCTION 

With the complementary investigations on heavy-ion col¬ 
lective flows and two precise measurements of heavy mass 
pulsars 0, the study of neutron stars (NSs) has become a 
more and more active field of research. Ab initio lattice QCD 
simulations and planned missions LOFT 0 and NICER will 
certainly promote further our current understanding of the un¬ 
derlying baryonic forces, high-density equation of state (EoS), 
and NSs’ core properties yj]. There might be three kinds of 
non-nucleonic components in the NS interior: free quarks 0 , 
mesons 0 , and hyperons fill]). But no ab initio calculations 
are available so far for their relevance and abundance in NS, 
because such calculations is unachievable due to the compli¬ 
cated nonlinear and nonperturbative nature of QCD. 

In the previous article 0 we have investigated the NS 
structure within the Brueckner-Hartree-Fock (BHF) approx¬ 
imation, which is currently one of the most advanced mi¬ 
croscopic approaches to the EoS of nuclear matter 0. In 
that paper BHF was combined with the confinement density- 
dependent mass model (CDDM) 0 for the quark phase to 
model hybrid stars (HSs), limiting ourselves to include only 
the confinement potential in the quark mass scaling. In this 
work we further extend our calculations by including also 
short-range leading-order perturbative interactions in the em¬ 
ployed quark matter EoS model (i.e., a new version of CDDM 
model 110]), and explore the consequences for HS structure. 
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Although in the high temperature region hadron-quark 
phase transition is crossover as explored by lattice QCD sim¬ 
ulations |ll], the order of phase transition at zero temperature 
is still an open problem, as well as the existence of a criti¬ 
cal end point (CEP) in the QCD phase diagram. We assume 
in the present work that the hadron-quark phase transition in 
cold NSs is a first-order one, and use the Gibbs construction 
to match the hadron EoS and the quark EoS for obtaining the 
mixed phase (see Refs. mm for more discussions on its 
dependencies on physical situations). That is, the pressure is 
taken to be the same in the hadron-quark mixed phase to en¬ 
sure mechanical stability, and would increase monotonically 
with baryon chemical potential. In the mean time, a global 
charge neutrality is assumed. Other authors used a smooth 
crossover |l4Ul5ll to obtain the transition. 

We provide a short overview of the theoretical framework 
and discussions of our results in Sect. II, before drawing con¬ 
clusions in Sect. III. 


II. FORMALISM AND DISCUSSION 
A. The hadron phase 

Let us first address the hadron phase, that is nuclear matter 
consisting of nucleons in [5 -equilibrium with electrons: 

n^p + e~ + V e - (1) 

We omit muons in the following calculations since they are 
irrelevant to the purpose of the present work. Under the con¬ 
dition of neutrino escape, this equilibrium can be expressed 
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as 

Pp = Pet (^) 

And the requirement of charge neutrality implies 

n p = n e - , (3) 

where /r, (n,) is the chemical potential (the number density) of 
component i. 

The chemical potentials of the non-interacting electrons are 
obtained by solving numerically the free Fermi gas model. 
The nucleonic chemical potentials required in Eq. ({2} are de¬ 
rived from the energy density of nuclear matter, based on the 
BHF nuclear many-body approach described elsewhere 0. 
Here the input bare nucleon force we employed is the Argonne 
V18 two-body interaction id, accompanied by a micro¬ 
scopic three body force constructed from the meson-exchange 
current approach [18]. The corresponding nuclear EoS repro¬ 
duces correctly the nuclear matter saturation point and fulfills 
several requirements from the nuclear phenomenology 
We mention here that the model developed in this article 
misses some important aspects, such as the inclusion of hy- 
perons. The interplay between hyperons and free quarks is 
quite important and deserves additional investigation, espe¬ 
cially when more reliable empirical inputs will be available, 
especially on the hyperon-nucleon and hyperon-hyperon in¬ 
teraction. 

Once the nuclear EoS and the nucleonic chemical potentials 
of nuclear matter are known, one can then proceed to calculate 
the composition of the hot ]3 -equilibrium matter by solving 
Eqs. © and ©, together with the conservation of the baryon 
number, n„ + n p = n g . Finally the total energy density £,y and 
the total pressure of the system are obtained after adding 
the standard contribution of electrons. 


B. The quark phase 

The quark phase is considered as a mixture of interacting u, 
d, s quarks in ]3-equilibrium with electrons: 

d ii T c T (4) 

s ^ u -\- e +V e , (5) 

s + u ^ d + u. (6) 

The crucial problem in studying the quark matter is to treat 
the quark confinement in a proper way. In the framework of 
the bag model, an extra constant, the famous bag constant B, 
is introduced which provides a negative pressure to confine 
quarks within a finite volume. That is, the quark mass is in¬ 
finitely large outside the bag, and finite and constant within 
the bag. As is well known, however, particle masses vary from 
the vacuum to a medium. Taking advantage of the density de¬ 
pendence, one can describe quark confinement without using 
the bag constant. Instead, the quark confinement is achieved 
by the density dependence of the quark masses derived from 
in-medium chiral condensates & That is the CDDM model 


we employed in the present study. A large amount of investi¬ 
gation have been performed in the framework of this model, 
and it has been developed greatly in recent years (see Ref. 0 
and references therein). 

In the employed CDDM model, strong interactions between 
quarks are mimicked by an equivalent mass to be determined: 

Hqcd : H k + 52 m qom + Hi (7) 

q=u,d,s 

= tik+ 52 

q—u,d,s 

where m q o (q = u.d.s) are the quark current mass, H\- is the 
kinetic term, H\ is the interacting part. The equivalent mass 
m q embodies all the interaction effects between quarks. That 
is, the contributions from both the scalar field and the Lorentz 
vector field can be included in this way (20f[. 

There are several ways to determine the equivalent mass in 
the literature (see [21] and references therein). Here we use a 
recently derived mass formula at zero temperature 0: 

m q = m q0 + mi = m q0 + + C/ij/ 3 . (8) 

n B 

where mi is the interacting mass, parameterized as a function 
of the baryon number density ng. D term is derived from the 
non-perturbative linear confinement of quarks (see also |9ll), 
and C term comes from the short-range leading contribution 
of perturbative interactions. These two terms correspond to 
the two leading terms in both directions when expanding the 
equivalent mass to a Laurant series of the holistic Fermi mo¬ 
mentum, respectively [HI]. The confinement interaction dom¬ 
inates at lower densities, while the perturbative interactions 
becomes more important at higher densities. This model gives 
reasonable results for the sound velocity. Also, strange quark 
matter in bulk still has the possibility of absolute stability for 
a wide range of parameters. As to the quark current mass, in 
our calculations we take m„o = o = 0 and m s o = 95 MeV. 

The parameter D has a lower bound I) 1 2 = 156 MeV, and 
an upper bound D 1 / 2 = 270 MeV l22h . The lower bound 
comes from the nuclear physics constraint, demanding that at 
P = 0, non-strange nuclear matter should be stable against de¬ 
cay to ( ud) quark matter. This leads to the condition E/A > 
M 56 p e c 2 /56 = 930 MeV for (ud) quark matter, which gives 
the above mentioned lower bound. The upper bound can be 
derived from a relation between D and the quark-condensate 
and the known range of values for this condensate lf22j] . The 
upper boundary of 270 MeV is in fact a very conservative 
one. According to the updated quark condensate determined 
nowadays very precisely by lattice QCD 0- a range of (161 
MeV, 195 MeV) can be obtained. Therefore in this work, 
we take two typical values of the confinement parameter as 
D 1 / 2 = 170 MeV, 190 MeV as inferred by the newest lattice 
QCD results ll23ll . 

The parameter C depends on how the strong coupling 
runs and it is determined so to have a_ upper bound of C = 
1.1676 0. Previous calculations 0 of pure quark stars 
employing a CDDM EoS lead to a maximum mass as high 
as 2 Mq with a parameter set (C,D 1//2 ) = (0.7,129 MeV). We 
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FIG. 1: (Color online) Extra chemical potential as a function of the 
baryon density, for both C = 0 and C = 0.7 cases with D 1 / 2 = 170 
MeV. 


FIG. 2: (Color online) Pressure of /3-stable quark matter as a function 
of baryon density for two values of C = 0,0.7, and two values of 
D 1 / 2 = 170 MeV, 190 MeV. 


will then employ C = 0.7 to perform calculations and change 
its value in a certain range for comparison. 

The relevant chemical potentials jx u , ji s , and ji e satisfy 
the weak-equilibrium condition (we again assume neutrino es¬ 
cape): 


Hu-lld = He, Hd = Hs- ( 9 ) 

The baryon number density and the charge density can be 
written as 

tiB = ^{>h, + n d +n s ), (10) 


q <2 = ^n u --n d --n s ~n e . (11) 

The charge neutrality condition requires qQ = 0. 

Since the quark masses m, are density dependent, the quark 
chemical potentials /q have an additional term iq with respect 
to the free Fermi gas model ( / = u.d.s): 


dsi d Vj 
d V, d iij 


+ £ 


de dnij 
dnij Bn, 


\J V, 2 + mj - jtq (12) 


The quark energy densities are 

e< = - bio ^/ p2+m ^ p2dp ^ (13) 

£ = ^£;, ( 14 ) 

i 


where v ( - = ( jt 2 /!/) 1 / 3 are the Fermi momenta and dnij/drn 
are derived from Eq. ® by taking the derivative of the baryon 
density. The quark pressure is calculated as p = —e + jv ft/'h- 
Solving Eqs. ©, © and (Ql}, the total energy density e q 
and pressure p q of the system can be obtained after adding the 
contribution of the leptons. 


In our previous calculations 0] there is only the term D 
in the quark mass scaling, which results in relatively low 
quark thresholds and small HS masses, similar to the calcu¬ 
lations |}24ll within the color dielectric model and the MIT bag 
model. As will be shown later, inclusion of the C term may 
bring strongly repulsive quark interactions, pushing the quark 
matter to appear at appropriately high densities, and conse¬ 
quently making massive HSs in the model. 

For a better understanding of this point, we present in Fig.Q] 
the modification of the extra chemical potentials jtq induced 
by the nonzero perturbative parameter C. The calculations are 
done with fixed Z) 1 / 2 = 170 MeV and two values of C = 0,0.7. 
We see that with increasing density, jtq always decreases. A 
nonzero C makes it decrease faster with the density. The de¬ 
crease is even more pronounced at higher densities, because 
the term Cn K is an increasing function of the density. Also, 
the extra chemical potential /q changes from positive to neg¬ 
ative values at high densities for C = 0.7, indicating that the 
term C brings repulsions, and quarks in this case are more 
strongly interacting with each other. This effect of C is oppo¬ 
site to that of D, since the latter arises from the quark confine¬ 
ment potential, and the increase of this parameter will bring 
attractions and soften the EoSs of the matter. Fig. [2] clearly 
demonstrates these effects, where the EoSs of /3-stable quark 
matter are shown for the two values of C = 0, 0.7, and the 
two values of D */ 2 = 170 MeV, 190 MeV. The increase of the 
confinement parameter D will soften the EoS, while the per¬ 
turbative parameter C will stiffen it. The repulsive nature of 
the term C will have crucial consequences for the structure of 
the resulting HSs as seen later. 


C. The mixed phase 

Let us consider the mixed phase made of nucleon matter in 
equilibrium with a gas of u, d,s quarks and electrons. Assum- 
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FIG. 3: (Color online) Quark fraction % as a function of the baryon 
density, for two values of C = 0,0.7, and two values of £)’/- = 170 
MeV, 190 MeV. 

ing the global charge conservation, the conservation laws can 
be imposed introducing the quark fraction % defined as 

X=V q /V, (15) 

where V is the total volume and V q the volume occupied by 
quarks. In terms of % the total baryon density, total electric 


charge and total energy density are written 


«B = (1 ~X) n N + X n q, 

(16) 

Qt = {l-x)QN + xQ q , 

(17) 

Et = (l-x)En + XE q , 

(18) 


respectively. The quantities «N(nq), (2n( Gq), and T’n (£q) are 
nucleonic (quark) number density, charge density, and energy 
density, respectively. 

Nucleonic chemical potentials are connected to quark 
chemical potentials as follows: 

AG = p u + 2p d , (19) 

= 2 Hu + Vd- (20) 

Therefore, there are only two independent chemical poten¬ 
tials. For a given total density «g, the two independent chem¬ 
ical potentials and the quark fraction % can be determined by 
solving the charge neutrality equation Qt = 0 and the pressure 
balance equation p n = p q - 

In Fig. [3] the quark fraction % is plotted as a function of the 
baryon density, for the two values of C = 0,0.7, and the two 
values of D 1 / 2 = 170 MeV, 190 MeV. As already known from 
our previous work, the larger D value pushes the threshold of 
quark matter phase to higher densities in the case of C = 0. 
This conclusion still holds when we adopt a nonzero value 




FIG. 4: (Color online) (Upper panel) Pressure of HS matter as a 
function of energy density and (lower panel) energy density and pres¬ 
sure of HS matter as a function of baryon density, for two values of 
C = 0,0.7, and two values of D 1 / 2 = 170 MeV, 190 MeV. The results 
of pure nuclear matter are also shown for comparison. 


of C in the quark mass scaling, increasing the critical density 
from 0.57 fm^ 3 to 0.64 fm^ 3 at (C,/) 1 / 2 ) = (0.7,190 MeV). 
Also, a nonzero positive C is found to have the same effect, 
when the C = 0.7 cases are compared with the corresponding 
C = 0 cases. A critical density 0.16 fm~ 3 (0.33 fm~ 3 ) with 
C = 0,T>'/ 2 = 170 MeV (190 MeV) is increased to 0.57 fm~ 3 
(0.64 fm~ 3 ) with C = 0 . 7 ,Z) 1 ^ = no MeV (190 MeV). 

In turn, when the C value is chosen to be negative, that is 
allowed by our model, quarks appear even earlier than the 
C = 0 case. This means that for typical model parameters 
free quarks could be present below the nuclear saturation den¬ 
sity 0.16 fm 3 , which is unphysical, so that we have to ex¬ 
clude negative C values in the present work. Actually, a pre¬ 
vious study showed that negative C values were allowed, for 
a simple inclusion of one-gluon-exchange interaction between 
quarks and degeneration 4 with respect to spin and flavor [|25]. 
This can be understood as follows: Since we start from the 
quark potential, the spin and flavor degrees of freedom should 









5 




r [km] 


mixed phase of the quark and nuclear matter. We then employ 
corresponding EoS models described above. They are shown 
in Fig. [4] where the pressures of HS matter as a function of en¬ 
ergy density are shown in the upper panel, and the energy den¬ 
sities and pressures of HS matter as a function of baryon den¬ 
sity are shown in the lower panel, for two values of C = 0,0.7, 
and two values of Z) 1 / 2 = 170 MeV, 190 MeV. These plots 
show again that the term C can stiffen the EoS, pushing the 
occurrence of free quark phase deeper into the star core, as 
discussed above. For the description of the NS’s crust, we 
have joined the hadronic EoSs above described with the ones 
by Negele and Vautherin [[299 in the medium-density regime 
(0.001 fm~ 3 < p < 0.08 fm~ 3 ), and the ones by Fe ynm an- 
Metropolis-Teller ll30ll and Baym-Pethick-Sutherland |3j3l for 
the outer crust (p < 0.001 fm 3 ). 

Fig- 0 shows the corresponding HSs’ masses as a func¬ 
tion of their central densities (upper panel) and also the HS 
mass-radius relations (lower panel) with the chosen param¬ 
eters. The results of the nucleon star are also shown for 
comparison. We see that a higher mass is achieved when 
the term C is included, for example, the value of 1.53M 0 
(1.65M 0 ) of the maximum mass obtained with C = 0,-D 1 / 2 = 
170 MeV (190 MeV) is increased to 2.1OM 0 (2.17M 0 ) us¬ 
ing C = 0.1,D 1 / 2 = 170 MeV (190 MeV). Only in the case 
of C = 0,D X ! 2 = 170 MeV, a pure quark core can be reached, 
while in other three cases, the most massive stars both have a 
mixed-phase core. We have also checked that, an even larger 
C value will bring heavier HSs, but always with a mixed-phase 
core. This means that no pure core is possible in the present 
HS model and quarks only appear in a limited region of the 
NS’s core. These results are consistent with a latest study us¬ 
ing the chiral effective field theory approach joined with the 
Polyakov—Nambu—Jona—Lasinio (PNJL) model |32jl. 


FIG. 5: (Color online) (Upper panel) HSs’ masses as a function of 
their central densities and (lower panel) HS mass-radius relations, 
for two values of C = 0,0.7, and two values of Z? 1 / 2 = 170 MeV, 
190 MeV. The results of the nucleon star are also shown for com¬ 
parison, as well as the recent 2-solar-mass constraint from the mass 
measurements of PSR J1614-2230 and PSR J0348+0432 d. 


be very important. Especially, the spin-spin interaction be¬ 
tween the quarks plays an important role in calculating the ef¬ 
fective repulsion [261- A more realistic inclusion of the quark 
potential as done in [27]] could then lead to a more reasonable 
result of the critical density around 0.5 fm 3 . Furthermore, 
when the term C is included, together with an upper quark 
threshold, the density range of the mixed hadron-quark phase 
is expanded (around twice with the chosen parameters). This 
should be a general result when the Gibbs construction is em¬ 
ployed to achieve a first-order phase transition. 

The stable configurations of a NS can be obtained from the 
well-known hydrostatic equilibrium equations of Tolman, Op- 
penheimer, and Volkov [ 28]. At variance with pure nuclear or 
quark stars, a HS may contain pure quark matter in the core, 
pure nuclear matter near the outer part, and, in between, a 


HI. CONCLUSIONS 

Summarizing, we have presented updated calculations of 
the transition from hadron to quark deconfined phase in NS 
matter, and also the HS structure based on our previous work. 
This extension concerns mainly the quark matter EoS, where 
we used a recent derivation of the quark mass scaling, includ¬ 
ing the leading-order perturbative interactions, in addition to 
the quark confinement. The derivation scheme allows us to 
modify largely the high-density behaviour of dense matter, re¬ 
sulting from a more repulsive quark interaction. 

We find that the quark thresholds are pushed to high densi¬ 
ties, together with large density jumps in the first order phase 
transition. Also, the EoSs are stiffened and the resulting HS 
maximum mass are shifted to higher values. Massive HSs as 
high as 2Mq are possible, consistently with two recent astro- 
physical observations of pulsars in binary systems. 

In the near future we plan to include the color superconduc¬ 
tivity since it is expected to play a role in dense quark matter at 
the density range discussed in the present work. Also, the ap¬ 
pearance of hyperons, missing in our present version of BHF 
model, could be studied in competition with free quarks, since 
a previous study lf33h . combining the same nucleon model 
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with the Dyson-Schwinger quark model shows that no hybrid 
star can exist if hyperons are introduced. Finally, we would 
like to study how such high-mass NSs are formed in a binary 
system. 
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